Coincidence measurement of residues and light particles in the reaction 56 Fe+p 

at 1 GeV per nucleon with SPALADIN 
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The spallation of 56 Fe in collisions with hydrogen at 1 A GeV has been studied in inverse kinemat- 
ics with the large-aperture setup SPALADIN at GSI. Coincidences of residues with low-center-of- 
mass kinetic energy light particles and fragments have been measured allowing the decomposition of 
the total reaction cross-section into the different possible de-excitation channels. Detailed informa- 
tion on the evolution of these de-excitation channels with excitation energy has also been obtained. 
The comparison of the data with predictions of several de-excitation models coupled to the INCL4 
intra-nuclear cascade model shows that only GEMINI can reasonably account for the bulk of col- 
lected results, indicating that in a light system with no compression and little angular momentum, 
multifragmentation might not be necessary to explain the data. 

PACS numbers: 25.40.Sc, 24.10.-i, 25.70.Pq 
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Spallation reactions play an important role in many 
domains ranging from astrophysics to intense neutron 
sources. Proton-induced reactions are also a way to study 
the de-excitation mechanism of a nucleus in a single hot 
source, and with less dynamical effects than in nucleus- 
nucleus collisions. They are often described as a 2-step 
model, with an intra-nuclear cascade (INC) phase fol- 
lowed by a de-excitation phase. Inclusive data on light 
particles emitted in the spallation process and, more re- 
cently, data on spallation residues, helped considerably 
in improving the models [1]. However, these are not suf- 
ficient to provide a real insight into the reaction mech- 
anism and the discrepancies observed between data and 
codes cannot be interpreted unambiguously with inclu- 
sive data. This is in particular due to the fact that the 
final observables are often both influenced by the cascade 
phase (especially by the remnant excitation energy) and 
by the de-excitation phase. A few more exclusive mea- 
surements exist but are generally limited to the study of 
the most violent collisions representing a small part of 
the total reaction cross-section (for a review see e.g. 0]). 

The need for a better understanding of spallation re- 
actions motivated the design of the SPALADIN setup 
at GSI, which aims at measuring in inverse kinematics 
and in coincidence all the spallation products with a low 
center-of-mass (cm.) kinetic energy, from neutrons to 
heavy residues. The restriction to low cm. energies, in 



fact due to geometrical acceptance limitations, largely 
favors the detection of particles from the de-excitation 
rather than the cascade phase. This allows to use the 
particle multiplicities as an indication of the excitation 
energy (E*) at the end of the cascade stage. 

The SPALADIN setup, partially described in @|, is 
based on the inverse kinematics technique where the ion 
beam is projected onto a liquid hydrogen target. The use 
of the large acceptance dipole magnet ALADIN permits 
to select the particles with a low cm. kinetic energy. 
Among other detectors, the setup comprises the large 
area neutron detector LAND, which provides neutron 
multiplicities, a time-of- flight wall and the multitrack 
and multiple-sampling time projection chamber (TPC) 
MUSIC-IV. The TPC enables a good charge identifica- 
tion for Z down to protons. The use of flash ADC's 
permits the simultaneous detection and reconstruction 
of several tracks. Its global efficiency was calculated to 
be 44 % for hydrogen, 78 % for He, 83 % for Li, and 
> 94 % for heavier fragments. The acceptance of the 
setup has been estimated through a GEANT 4 simula- 
tion as being around 20 % for protons, 80 % for helium 
and > 97 % for heavier fragments. Our relatively poor 
efficiency and acceptance for protons is the reason why 
data on protons are not shown here. Data are obtained 
after empty target subtraction, once normalized to the 
number of incident beam particles counted with a thin 
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plastic scintillator also used as a trigger. In the following, 
they will be compared to the codes filtered by GEANT 4. 

The first experiment using this setup was performed on 
the system Fe+p at I GeV per nucleon. In addition to 
the coincidence data, inclusive element cross sections and 
recoil velocity distributions were obtained, which will be 
presented in more details in a forthcoming paper. The 
element cross sections obtained with SPALADIN agree 
with previous data from the FRagment Separator (FRS) 
0, El (not shown here) except for carbon and nitrogen, 
which deviate by 25%. 

As was done for the FRS data comparisons with 
several de-excitation models associated with the intra- 
nuclear cascade INCL4 [6] allows to draw some first con- 
clusions (see Fig. [1]). The combination of INCL4 with 
the de-excitation model ABLA Q, which is generally 
rather successful in predicting data in heavy systems, se- 
riously underpredicts the production of fragments with 
charges between 3 and 10 (hereafter denoted as Interme- 
diate Mass Fragments or IMF). This can be ascribed to 
the fact that ABLA is a "classical" de-excitation model 
in which only the evaporation of neutrons, protons and 
alpha particles is possible. To account for the emission 
of IMF, several de-excitation mechanisms can be envis- 
aged: A generalized evaporation, as modeled in GEM Q 
which evaporates particles up to Mg; a binary splitting 
as described in GEMINI M, which uses the Moretto-type 
transition state model [1 01 ] for the emission of intermedi- 
ate mass fragments down to a certain limit, the evapora- 
tion of lighter fragments being treated within the Hauser- 
Feshbach formalism; the ope ning of multifragmentation 
channels, as included in SMM 11]. Default parameters of 
the models have been used except for GEMINI for which 
the transition from the Hauser-Feshbach to the Transi- 
tion State Model is made for Z>4 as recommended by 
the author We can see that IMF cross-sections are 
largely under-predicted by GEM while GEMINI better 
agree with the data and SMM overestimates them. Sim- 
ilar conclusions were drawn also in 13j. However, from 
these inclusive data it is not possible to draw more precise 
conclusions on the actual mechanism. 

We cannot exclude that the underprediction of light 
elements could be due to the intra-nuclear cascade model, 
INCL4, which, for instance, could predict a too small 
excitation energy of the de-exciting nucleus. However, it 
was shown in [4| that a cascade model leading to higher 
excitation energies would certainly help populating the 
IMF region but at the detriment of the high Z part of 
the spectrum. It was also noticed that excitation energies 
given by INCL4 at the end of the cascade are remarkably 
close to those given by the ISABEL code [1J] , which uses 
a very different modeling of the cascade. 

With SPALADIN, it is possible to decompose the ele- 
ment production cross section into the different reaction 
channels. The decomposition has been done according 
to the numbers of Z=2 particles and fragments (Z>3) in 
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FIG. 1: (Color on line) Element production cross sections 
compared to ABLA, GEMINI, SMM and GEM, associated 
with INCL4. 



each event. The results are presented in Fig. [2] (upper left 
panel) . It can be seen that the heaviest fragments (down 
to charge 18) are most often detected alone (no Z=2) (in 
fact they are accompanied by only neutrons or hydrogen 
isotopes) while lighter ones are generally associated with 
helium emission. As expected, production cross sections 
of IMF are mostly populated by events with at least two 
fragments. This explains why in Fig. [T] this part of the 
spectrum is badly reproduced by classical evaporation 
codes as ABLA. The same decomposition has been done 
with the de-excitation codes mentioned above, all associ- 
ated with INCL4 (three other panels of Fig. ^ . They are 
compared to the experimental data represented by the 
solid lines. As seen earlier, the generalized evaporation 
model, GEM, predicts too little IMF production. The 
interpretation for production cross sections of fragments 
with a charge between 10 and 15 is more delicate. In the 
experimental data, this part of the spectrum is mostly 
dominated by events with at least one helium particle in 
the final state. The fact that GEM underestimates these 
cross sections could be due to a lack of helium production 
at high E*, as will be shown in the following. 

The case of SMM illustrates that the data obtained 
with SPALADIN are much more constraining than in- 
clusive production spectra: For instance, in the region 
8<Z<15, the total production cross sections are rela- 
tively well reproduced by SMM. Actually, this results 
from an overprediction of events with 3 fragments com- 
pensated to some extent by an underestimation of events 
with helium emission. Finally, we can see that the best 
agreement is obtained with GEMINI that reproduces 
rather well the different contributions although events 
without Z=2 are a little overestimated. 

In order to understand the transition from classical 
evaporation to IMF emission, we have looked for vari- 
ables that could be related to the excitation energy at 
the end of the cascade stage, E* . Two different and 
nearly complementary variables have been chosen: The 
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FIG. 2: (Color online) Contributions of the different channels 
to the element cross section, measured (upper-left panel) and 
calculated with GEM, SMM and GEMINI, associated with 
INCL4 and filtered through GEANT4. The solid lines delim- 
iting the different regions are taken from the upper-left panel 
and shown on the other panels. Fragment means Z>3. 



multiplicity of neutrons plus Z=2 particles, as we mostly 
detect de-excitation particles (and poorly Z=l), and the 
variable Zb oune i. The latter is defined as the sum of all 
charges strictly larger than one, as originally introduced 
by the ALADIN collaboration [15[ > and is related to the 
charge remaining bound into fragments. The choice of 
one or the other variable is dictated by the necessity to 
avoid possible correlation effects when looking at certain 
observables as a function of one variable. Their correla- 
tion with E* has been studied with the different models 
filtered by GEANT4 and found to be practically inde- 
pendent of the de-excitation model. The case of Zhound 
can be seen on Fig. [3] left. This indicates that the range 
of E* associated with a given bin of Zhound, for instance, 
can be determined rather unambiguously. Of course, the 
distribution of E* over all events could depend on the 
INC stage. The distribution of E* divided by the mass 
of the remnant An is displayed in Fig. [3] right. Events 
with E*/Ar > 4 MeV (often claimed to be the threshold 
for multifragmentation 0]) represents 6.5% of the cross- 
section. Events with at least 2 fragments are also shown: 
with SMM (dotted curve), as expected for a multifrag- 
mentation model, they are associated with the highest 
E*/A R , beginning around 3 MeV, while with GEMINI 
(dashed curve) they are more broadly distributed. 

The mean multiplicities per event of the different types 
of fragments are shown as a function of Zhound on Fig. [4] 
It reveals that mean multiplicities of Z > 2 fragments 
remain low, even at high E* . Actually, events with 2 
(resp. 3) fragments represent 32±3 (resp. 2±0.2) mb 
out of a total cross section of 777±79 mb. With GEMINI 
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FIG. 3: (Color online) Left: E* at the end of the cascade 
versus the variable Zbound simulated for the different de- 
excitation models, after GEANT4 filtering. The vertical bars 
indicate the standard deviation of the distributions in the case 
of GEMINI. Right: Distribution of E* /Ar for all events (full 
line) and for events with at least 2 fragments for GEMINI 
(dashed line) and SMM (dotted line) . 
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FIG. 4: (Color online) Mean multiplicities as a function of 
Zbound (explanations in text). Lines for the codes are stopped 
when the number of events for a given Zbound becomes lower 
than 30. 



(resp. SMM) the cross-section of events with at least 2 
fragments is 50 (resp. 80) mb. An apparent odd-even 
effect is observed in the mean multiplicity of Z = 3 and 
Z = 5. This is linked to the mainly binary character of 
the breakup and to the fact that it is unlikely to produce 
two fragments with both an odd number of protons. 

The generalized evaporation of both GEM and GEM- 
INI, respectively through the Weisskopf and Hauser- 
Feshbach formalism, reproduces very well the mean mul- 
tiplicities for Z = 3. For Z = 4, GEM is still very good 
but GEMINI less. For heavier fragments, however, the 
multiplicities predicted by GEM remain much too low, 
whereas they are well reproduced by GEMINI. The com- 
parison with SMM shows that it globally strongly over- 
predicts the IMF emission. As shown in [J], the contribu- 
tion of multifragmentation in SMM is about 30% in the 
region Z < 10, however the use of a breakup volume equal 
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to 3 times the pre-fragment volume (code default value) 
is questionable in the case of spallation. In fact, the IMF 
multiplicities remain globally the same when the breakup 
volume is reduced to the pre-fragment volume. In that 
case, the lower number of IMF emitted through multi- 
fragmentation is globally compensated by an increasing 
number of IMF emitted through the SMM generalized 
evaporation module. The case of helium emission is also 
very interesting because it keeps increasing when the ex- 
citation energy increases. This is well predicted by SMM 
and GEMINI but not by GEM, which predicts a plateau 
for Zb oun d < 20. This behavior is linked to the compe- 
tition between neutrons, protons, He and IMF emission 
at high excitation energy, which depends on the model- 
ing of the Coulomb barrier, the density of states or the 
transmission coefficients in the different codes. 

In order to better understand the mechanism responsi- 
ble for the IMF emission, we have looked at the difference 
between the highest charge, Zl, and the second high- 
est charge, Z2, for multi-fragment events (Fig. [5j . The 
events were divided into 3 bins in particle multiplicity, 
going from low to high excitation energies. A transition 
from an asymmetric breakup to a more symmetric one 
is observed when E* is increasing. While GEM rather 
well reproduces the case of large Zl — Z2, it is totally 
unable to predict more symmetric breakup channels, as 
expected from a model with only evaporation channels. 
The evolution with E* is better reproduced by SMM but, 
as mentioned before, the global IMF emission is too high 
and the shape at low E* not correct. This could be due to 
the evaporation part of SMM. Only GEMINI reproduces 
very well the results, in shape and cross-sections. These 
observations indicate that evaporation, even generalized 
to IMF emission cannot account for the observed sym- 
metric breakups and that a formalism like the Transition 
State Model could explain the production of fragments 
with Z > 4 better than a model with multifragmentation, 
even for events corresponding to E* /Ar ~ 4 MeV. This 
does not support the conclusion of [H| based on veloc- 
ity measurements that IMF are formed in simultaneous 
breakup decays. 

To conclude, the simultaneous measurement of light 
particles and spallation residues, performed in inverse 
kinematics with the SPALADIN setup, on the system 
Fe+p at 1 GeV per nucleon, has allowed for the first 
time the decomposition of the total reaction cross-section 
into the different decay channels and the study of its 
evolution with excitation energy. The comparison with 
the predictions of different de-excitation models coupled 
to INCL4 confirmed that the observed high production 
cross-sections of IMF, already observed in inclusive mea- 
surements, cannot be explained by standard evapora- 
tion models, even generalized to the emission of IMF, 
as GEM. This model does not predict the appearance 
of symmetric breakups observed when the excitation en- 
ergy increases. The use of models with other mecha- 
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FIG. 5: (Color online) Charge difference between the 2 heavi- 
est fragments detected in multi-fragment events, for 3 bins in 
the multiplicity of neutrons and Z=2 particles, corresponding 
to average E* /An of respectively 3.1, 3.8 and 4.5 MeV. 



nisms for the IMF production, such as multifragmenta- 
tion in SMM or the Transition State Model asymmetric 
splitting modes of GEMINI, helps to reproduce the data. 
It appears that the best description of the whole set of 
coincidence data is obtained with the INCL4+ GEMINI 
code, in particular as regards the dependence with exci- 
tation energy. It seems therefore that in a light system 
as Fe+p, with no compression and little angular momen- 
tum (9h in average) multifragmentation might not be 
necessary to explain the data. The better understanding 
of the reaction mechanism reached through coincidence 
measurements of the present kind will certainly allow the 
development of more reliable simulation tools for appli- 
cations of spallation reactions. 
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